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ABSTRACT 


A  six  degree  of  freedom  model  for  di-  >r  di¬ 
gital  simulation  of  the  trajectory  of  an  -  an 

unguided,  fin-stabilized  rocket  is  developed.  sveloped. 

A  derivation  of  the  equations  and  an  expla-  expla¬ 
nation  of  the  coordinate  systems  are  pre-  pre¬ 
sented.  The  development  assuses  that  the  .  the 

trajectory  will  be  over  a  rotating  planet  anet 

with  a  variable  ataosphore.  A  space- 

variable,  three-dimensional  wind  vector  is  .or  is 

assumed. 

The  equations  of  motion  are  derived  ved 

from  Newton's  Laws  of  Motion.  The  aero-  :ro- 

dynamic  forces  and  moments  are  based  on  the  on  the 

theory  of  stability  derivatives  and  the  he 

assumption  of  linear  aerodynamics.  The  lie 

body  axes  are  assuned  tr  be  principal  axes  axes 

of  inertia. 
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INTRODUCTION 


Various  unguided  rocket  systems  have  been  used  at  White  Sands 
Missile  Range  (WSMR)  for  several  years  for  high-altitude  research. 

Present  planning  indicates  that  an  increased  number  and  variety  of  such 
rockets  will  be  fired  at  WSMR  for  novel  experimentations  and  for 
re-entry'  studies.  The  use  of  the  unguided  rocket  is  desirable  since 
it  is  usually  less  complex,  and,  therefore,  more  reliable;  and  is  also 
less  expensive  than  one  with  guidance.  The  principal  disadvantage  is 
the  inherent  dispersion,  since  there  is  no  guidance.  The  largest 
contribution  to  the  dispersion  of  unguided  rockets  is  due  to  the  wind 
encountered  during  flight.  This  dispersion  can  be  made  less  signifi¬ 
cant  by  applying  a  procedure  for  determining  the  wind  effect  and  then 
adjusting  the  launcher  so  as  to  compensate  for  the  wind  effect,  there¬ 
by  achieving  the  desired  trajectory. 

A  real-time  prelaunch  impact  prediction  system  has  been  developed  4 

for  use  at  WSMR  to  compensate  for  the  effect  of  the  wind  [1,2].  An 
integral  part  of  this  system  is  the  trajectory  simulation  equations. 

The  equations  presently  used  by  the  system  were  developed  by  Dr.  Everett 
L.  Walter  [3].  This  report  presents  the  development  of  a  new  set  of 
equations  which  will  be  evaluated  for  possible  replacement  of  the 
Walter  equations. 

The  development  of  a  trajectory  simulation  model  rin  fellow  numer¬ 
ous  approaches  and  ctn  vary  from  quite  simple  to  extremely  det^.^ed, 
depending  primarily  upon  the  purpose  for  which  the  model  is  to  be  used. 

Since  this  model  is  to  be  used  solely  for  flight  simulation  in  a  real¬ 
time  system  to  determine  wind  effect  on  rockets  having  range  and  alti¬ 
tude  of  less  than  600  miles,  the  following  assumptions  are  made: 

a.  The  rocket  is  assumed  to  be  a  rigid  bodf  with  six  degrees 
of  freedom. 

b.  Tha  body  axes  are  principal  axes  of  inertia. 

c.  Linear  aerodynamics  are  adequate  for  determining  the  aero¬ 
dynamic  forces  and  moments. 

d.  The  earth  is  a  sphere. 

e.  Gravity  follows  the  inverse- square  law. 

f.  The  thrust  vector  acts  parallel  to  the  longitudinal  axis 
of  the  rocket. 
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LIST  Ul  SYMIS'TS 


X’.Y'.Z* 

»*> 

*»y  t* 

TX2x 
^X*^Y * 

kx- *^Y* **z' 

k  ,k  ,k 
x  y  z 

It 

<•*1 

P.q.r 

w 


♦l 

♦l0 

*1 »12»13 


"l  *m2*ia3 


•l*  2*"3 


Cl 


Launcher  coordinate  system 
Inertial  coordinate  system 
Body  coordinate  system 
Linear  transformation 

Unit  vectors  along  the  axis  denoted  by  the  subscript 

Position  vector  of  missile  center  of  gravity 
in  inertial  system 

Rotation  of  x,y,z  system  in  X,Y,Z  system 

x,y. z  components  of 

Earth's  rotational  velocity 

Angle  determined  by  the  projection  of  R  in  the 
X,Y  plane  and  the  X-axis 

Latitude  after  launch 

Initial  Latitude 


Direction  cosines  of  the  x,y,z  axes  respectively 
with  respect  to  the  X,Y,Z  coordinate  system 


Aerodynamic  stability  derivatives 


3C, 


'N 


3  a 


where 

C^>  is  the  normal  force  coefficient; 


U.V.W 

ZF  ,ZF  ,ZF 
x  y  z 


MaSj  of  rocket 

Atmospheric  pressure 

Area  of  exit  nozzle  of  rocket  motor 

Thrust  as  measured  by  a  static  test 

Atmospheric  pressure  at  static  test  site 

< 

Acceleration  due  to  gravity 

Value  of  g  at  sea  level 

Mean  sea  level  radius  of  earth 

Height  above  mean  sea  level 

Elevation  angle  at  launch 

Azimuth  angle  at  launch 

Total  force  acting  on  the  uissilo 

Sum  of  the  external  moments  acting  on  the  missile 

X,Y,Z  components  of  Mt 

Missile  velocity  vector 

Angular  momentum  vector 

V  referred  to  the  inertial  system 

V  referred  to  the  body  system 
x,y,z  components  of 

x,y,z  components  of  F 
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i 


x,y,~  components  of  .1 


I  -I  -I 


XX 

xy 

XZ 

-1 

xy 

i  -i 

yy  y: 

Tensor  of  inertia 

-I 

-i 

i 

x  7  y 


. iv  .Vim!  cm  portents  in  the  launcher  system 

x* 1  v  :•  1 

;;  .W  ,K  i ml  components  in  the  hotly  system 

V  V  7 


V  Speed  of  missile  with  respect  to  the  wind  vector 

a 

u* »v*,w'  X,y,z  components  of  $a 

Angle  of  3ttack 
Angle  of  sideslip 
Absolute  .ingle  of  attack 
Auxiliary  angle  of  attack 
1  tary  angle  of  s  ides  lip 
Speed  of  sound 
Atmospheric  density 
Dynamic  pressure 
Reference  area  for  aerodynamics 
d  Reference  length  for  aerodynamics 

COORD!  X  ATI;  SYS  TKMS  ANP  TRANS  FORMATIONS 

Three  right-hand  orthogonal  ccordinatc  systems  are  used  to  describe 
the  missile's  position  i si  space.  They  are  the  launcher  coot uinntc  sys¬ 
tem  with  X*  axes,  the  missile  coordinate  system  with  x  ,y  ,z  ;:;;ss 

and  the  inertial  coordinate  system  with  X,Y,Z  axes. 

The  X' ,Y',Z'  system  has  its  origin  at  the  launcher  and  rotates  with 
the  earth.  The  positive  X'  axis  points  cast,  the  positive  Y'  axis  points 
north,  and  the  positive  Z*  axis  points  outward  along  the  radius  vector 
from  the  center  of  the  earth. 


3 

<3 


Ihe  x,/,;  systci.  is  ;i  i.Mvliiy  sy^ivr:  i.iti:  its  crigut  at  the  center 
of  gravity.  Cg,  ot  the  Missile  il»e  x  axis  coincides  with  the  longi¬ 
tudinal  axis  of  tlie  rochet  am!  is  positive  toward  the  nose  Let  6  he 
the  angle  between  the  x  axis  ami  the  Z'  axis  Measure!  iron  the  Z’  axis. 
The  y  axis  lies  in  the  X’,Y‘  plane  ami  is  positive  in  the  direction  of 
positive  8 .  The  positive  z  axis  is  chosen  so  that  we  have  an  orthogonal 
right-hand  system. 


Tlxe  X,Y,Z  system  lias  its  origin  at  die  center  of  the  earth.  The 
system  is  oriented  so  that  the  X  and  Y  axes  lie  in  the  equatorial  plane 
with  the  Y  axis  initially  passing  through  the  longitude  of  tne  launcher. 
The  Z  axis  is  coincident  with  the  earth's  axis  and  positive  toward  the 
North  I’ole.  This  system  does  not  rotate  with  the  earth. 

The  linear  transformations  between  these  systems  will  be  denoted 
by  e.g.,TX2x  where  the  left  hand  subscript  denotes  the  domain  of  the 
mapping.  To  derive  those  transformations,  let  kv,ky ,k2,kxi  ,kv*  •  *2' » 
kx.k  kj  be  unit  vectors  along  the  axis  denoted  by  the  subscript.  Let 
R  be  position  vector  of  the  C„  and  lot  define  the  rotation  of  the 
x,y,z  system  in  the  X,YJ,Z' system  with  x,y,z  components  denoted  by  p,q, 
and  r  respectively.  Then 

R  ■  Xkx  ♦  Yky  ♦  3fcz*  (1) 


and 


“M  •  Pkx  *  t|ky  *  rV 


'2) 


It  can  be  seen  from  Figure  1  that  the  direction  cosines  of  R  in  the 
X, Y , Z  system  are  (cost|  sinv,  cos*L  cosy,  sinfL)  where 

sin  y  *  X(X2  ♦  Y2)  ’/2, 

cos  y  ■  Y(X2  ♦  Y2)*^2, 

sin4L  =  Z(X2  ♦  Y2  ♦  Z2r>/2,  am! 

co$4^  *  (X2  ♦  Y2)1^2  (X2  ♦  Y2  ♦  Z2)*1/2.  (3) 

Also,  from  Figure  1,  the  transformation  Tx2X*  ls  h*ven  by  the 
following  equations: 


kX'  *  *cos(wAt)kx  - sin(uAt)k^ 

k  (  ■  sin$  sin(uiit)k  ■  sine  cos(u£t)k  ♦  cose  k 
Y'  Lo  X  Lo  Y  t>  1 

k_,  *  -cos$  sinfwutlk..  *  cos?,  cos(wat)ky  ♦  sin$.  k_.  (4) 
"  o  A  o  **o  “ 


0 


l'he  inverse  •  "*x  t'lC  trans formation  ,y,  is  given  by: 

ky  -  * c us ( um t )  k  y ,  ♦  sin'|^  sin(i.ut)kyi  -  cos^i^  sin(w&t)k~, 

ky  *  -sintiuitjkj.,  -  sin<f Lq  cos(uAt)ky,  ♦  cos^  cos(uAt)k2, 


k,- 


■  cos 


*1. 


0 


ky  *  *  5111$^  k,,. 


(5) 


(observe  that  since  cadi  of  the  coordinate  systems  is  orthogonal  the 
inverse  of  3  transfornation  is  just  the  transpose.) 


Let  (lj  «lj *13) 
tion  cosines  or  the 
transformation  TX2x 


,  (r.,  ,m2,ra3)  and  (nItn,»aj)  be  the  respective  dircc- 
x,  y,  and  z  axes  in  the  X,Y,Z  system.  Then  the 
and  its  inverse  Tx2x  arc  given  by  the  sets  of  equations 


and 


kx  *  XlkX  *  *2kY  *  l3kZ 
ky  '  mtkX  *  B2kY  *  n3kZ 
kZ  "  "lkX  *  n2kY  *  n3kZ  » 


kX  * 


kY  * 


ljkx  ♦  nijky  *  n,k2 
l2kx  *  nipky  *  i^i- 

l  jkjj  *  n3ky  ♦  njkj  | 


respectively. 


(f) 


(7) 


Differential  equations  for  ’  i,  mi,  11  j,  i  ■  1,  2,  3  are  obtained 
as  follows:  let  the  linear  transfornation  be  expressed  by  the 
matrix 


let  the  matrix  of  the  linear  transfornation  be  expressed  by 

(see  Lquition  (4)): 
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*  COSuiU* 


-sinuit 


>'  =  slnCl.0si*‘<-it  -sinc^cosuit  c°s^|  q  I  ; 
\  ■C0S$L0sin“‘it  cos^  cosw&t  si n$j/ 

the  transformation  Tx2}.  is  expressed  by  the  matrix 
/  h  ml  ni  \ 


(h  m,  nj  \ 

l 2  m2  n2  I  . 

*3  n3  n3  / 


Hence  A  -  DC.  Thus, 


•  •  «  «  .  • 

A  -  DC  ♦  DC  ■  L'U  A  *  DC  . 


Now  to  determine  the  elements  of  A  it  suffices  to  determine  the 
elements  of  C  and  express  these  in  terms  of  A.  Consider  a  rotation, 
d*.  of  the  x,y,z  system  in  the  X'.Y^Z*  system.  The  unit  vector  k 
changes  by  an  i,nount  of  dkx.  The  projection  of  dkx  in  the  x,y  plane 
is  ky •  dkx  (see  Figure  2). 


ky-dkx 


FIGUR1:  2 

For  small  rotations,  i.e.,  for  small  changes  in  time,  the  components 
of  rotation  about  tho  z  axis  are  given  by 

d*x  ■  sindi^z  •  ky'dhx. 


r  ■  dijij/ dt  ■  ky*dkx/dt  ■  ky'kx. 


i 


Similarly 


p  *  kz  *  k.y 

ami 


(10) 


4 


kx*kz 


(11) 


Now  kx  lies  in  the  y-z  plane,  ky  in  tlie  x-t 
x-y  plane.  Thus  1 


plane,  and  kz  in  the 


and 


kx  3  kyA3  ”  kz**7» 


Z  2‘ 


“  kjAj  -  kxBj, 


(12) 


kz  “  kxB:  ’  ky*V 
Since  the  system  is  right  handed 


Hence, 


kx  ■  ky  *  kz* 
ky  “  kz  K  kx* 
kz  "  kx  *  V 


ky  "  kZ  *  ky  *  kj 


or 


(13) 


(14) 


ky^3  "  kz^2  "  (kz^l  “  kx®3^  *  kz  *  ky  *  (kx®2  "  ky®i) 

-  (-kx  X  kj)B3  *  (ky  x  kx)Bz  -  kyBj  -  kzB2.  (IS) 

So  A3  ■  B3  and  A,  ■  U2.  Similarly  At  ■  Bj.  From  liquations  (9),  (10), 
(12),  and  (13),  we  have 

A,  -  P,  A2  -  q,  A,  »  r.  (10) 

Now 

k,  «  i,kx,  ♦  i2ky,  ♦  i3kz, 

S'  *  ^iSc'  *  ’■'zS'  *  hkz'  (17) 

•  •  •  ■ 

kz  *  clkx'  *  c2ky'  *  c3k,  t  • 


« 


0 


"iiius  i'r>M.-  i.quatioiis  to) , 

(12).  (bl,  (10),  aikl  (17) 

v.c  have 

"•Y  +  “zV  * 

;ljk..  =  k,  =  Irbj-qcJk  1 

*  (rh2-qc2)ky. 

♦  (rb3-qc3)k. 

f  • 

So  a  =  raj  -  qb^ 

i  ■  1.  2,  3. 

(18) 

Similarly 

hi  -  pci  - 

i  *  1.  2,  3, 

(ia) 

and  ci  *  qa^  -  pb3 

i  -  1.  2,  3. 

Now  from  the  matrix  equation  A  -  DC  we  can  determine  a*,  6^  and 
i  »  1,  2,  3  in  terms  of  13,  ■£,  and  n^  i  *  1,  2,  3  Performing  this 
substitution  and  the  matrix  algebra  and  differentiation  operations  in 
equation  (3)  we  obtain,  after  simplification, 

1|  *  rail  •  qnj  -  ulj 

*  rm2  •  qn2  +  wlj 

13  *  r»3  -  qn3 

ml  *  Pnl  *  til  -  “®2 

m.  »  pn2 

“3  *  Pn3 

nl  *  *1^1  “  P*j  "  wn2 

n;  *  ql2  *  pm 2  *  “"l 

n3  -  ql3  -  pOj  . 

INITIAL  CONDITIONS  POP  TIIL  x,y  ,l  SYSTEM 

Initially  the  x  axis  and  the  Z*  axis  font  an  angle  60  (this  is  the 
elevation  of  the  x  axis  angle  and  is  measured  from  Z'};  and  the  projection 
of  the  x  axis  in  the  X'Y'  plane  forms  an  angle  oq  with  the  Y*  axis  (this  is 
the  azimuth  angle  of  the  x  axis  and  is  measured  clockwise  from  ”’).  Since 
the  y  axis  lies  initially  in  the  X*  ,Y’  plane,  it  fonts  an  angle  Oq  ♦  ‘JO  with 
the  Y'  axis  (see  Figure  3).  Thus,  initially, 

kx  «  kx,  sin«0sina0  ♦  ky,  sin0ocosaQ  ♦  kz«  cos0Q, 

ky  ■  kxi  cosao  ♦  ky»  sinaiQ, 

k2  *  kx,  sinaoCosSo  ♦  ky  cosaucos60  -  kz  siite0,  (21) 


rlj  ♦ 


rl, 


(20) 


10 


Fran  Equation  (5)  we  have,  initially, 

kx  r  .kx,  , 

ky  -  -sinfLo  kY,  ♦  coa*Lo  kj, 

kz  «  cos$L  kY>  ♦  sin*^  kz»  .  (22) 


Tlius  the  initial  conditions  for  l^,m^,n^,  i  *  1,2,3,  are  j,ivcn  by 


r 

M 

CM 

«-e 

sin0osinao  sin0ocosao  cos60 

”-l  0  0 

iHj  n2  mj 

■ 

cosaQ  ’*ina0  0 

0  -sin$,  cos*. 

Lo  Lo 

nl  n2  n3 

sir.aocos0o  cosaQcos90  -sir.90 

0  cos*L<>  sin4L(j 

_ 

—  — 

-  — 

DERIVATION  OF  HIE  EQUATIONS  OF  MOTION 

n»is  section  presents  the  derivation  of  the  equations  of  notion  of 
a  rigid  body  in  inertial  space.  This  derivation  can  be  found  in  several 
references  [4, 5,0,7!  and  is  presented  here  for  completeness  only.  The 
two  basic  equations  which  define  the  motion  are,  by  Newton's  1  aws. 


11 


I 


3£  «,iV1 


•'t  •  Jt  (il) 


These  equations  state  that  the  sum  of  the  applied  forces,  P,  is 
equal  to  the  time  rate  of  change  of  the  linear  momentum,  m$,  and  the  sum 
of  the  external  moment  Mt,  is  equal  to  the  time  rate  of  chance  of  the 
angular  momentum,  (t 

HIE  TRANSLATIONAL  ACCELERATIONS 

The  differential  equations  describing  the  motion  of  a  missile  sys¬ 
tem  (Equations  (24)  and  (2SJ)  must  always  be  referred  to  the  inertial 
system  (the  X,Y,2  system)  for  solution.  This  is  not.  meant  to  imply  that 
every  computation  is  done  in  the  inertial  system--* o  the  contrary  most 
computations  are  performed  elsewhere.  It  does  mean  that  the  basic 
equations  of  motion  must  be  referred  to  the  inertial  system  and  compu¬ 
tations  in  other  frames  related  by  appropriate  transformations. 

Let  V1  be  the  missile  velocity  in  the  inertial  system  and  the 
velocity  in  the  body  system  (throughout  this  section  the  superscripts 
I  and  b  will  be  used  as  above).  It  is  clear  that  V1  ■  V*>.  Since  Newton's 
laws  of  motion  are  valid  only  in  an  inertial  system,  the  derivatives  of 
(22)  and  (23)  must  be  in  the  inertial  system.  Now 

iir  c?I)  ■  nr  (**)  .  r*: 

The  right-hand  side  of  (26)  is  mixed  and  cannot  be  worked  >ith  as 
it  stands.  However,  the  derivative  with  respect  to  time  in  the  inertial 
system  can  be  written  as  the  operator 


Ut  "  Jt  *  V 


Hence , 


r,  <*“>  ■(£  *  -,■>*  ■  £  •  v 


Expanding  the  above  differential  equation  gives 


>  .  d 


_  (uhx  *  vk  ♦  wkz) 


ky  kz 


cr 


~*u  ,  , 

V  =  (u  ♦  qw  -  rv)kx  ♦  (v  ♦  ro  -  pw)ky  ♦  (w  *  pv  -  qu)kr»  (30) 
Now  from  the  force  equation, 

EFb  *  ,  (31) 

wc  have  translational  equations  of  motion: 

EPX  *  m(u  ♦  qw  -  rv)  , 

ZFy  ■  m(v  ♦  ru  -  pw) , 

IPj  >•  m(w  ♦  pv  -  qu)  ,  (32) 

where  £FX,  £Fy,  and  IP2  are  the  applied  forces  in  the  body  coordinate 
system. 

Tllli  ROTATIONAL  ACCCLL RATIO]  S 


The  relationships  expressing  the  rotational  notion  are  obtained  in 
a  straightforward  manner.  The  components  considered  come  from  two  sources, 
the  time  rate  of  change  of  the  moment  of  momentum  and  the  external  iy 
applied  moments.  Hie  moment  of  a  rigid  body  about  its  center  of  mass  is 
given  by: 


*xx 

-Ny 

-i  1 
xz 

11 

y 

at 

_Ixy 

lyy 

’lyz 

-»z  . 

-I 

U  V  -7 

-i 

v  y 

lzz  - 

P 

q 


(33) 


Since  the  body  coordinate  axes  are  assumed  to  be  principal  axes  of 
inertia,  the  off  diagonal  terms  of  the  3  by  3  n.-trix  arc  equal  r.o  zero. 
Thus  we  can  rewrite  liquation  (33)  as 

n  *  Ixxpkx  ♦  Iyyqky  ♦  Izzrkz  .  (34) 

Thus  we  obtain 

11  ■  ‘ScUxxP  +  1XXP  ♦  <xzz  ■  ^yjq*-]  *  ky[Iyyq  *  lYYq  * 

dxx  '  lyy)Prl  *  M*zzr  *  liz'T  *  dyy  -  IxxJPql '  (35> 

where  the  differentiation  of  moment  of  inertia  refers  to  that  change  at 
constant  mass  only.  The  change  in  angular  momenti*  due  to  the  escaping 


gases  :  l  e  .  Juc  to  change  m  muss!  is  considered  in  the  external  moments 
as  th  jet  damp mg  tionent, 

Now,  if  we  denote  the  sums  of  the  x.y.z  components  of  the-  external 


moments  by  L  >1,  am 

i  (to 

adopt 

stamia  n) 

L  =  lxxV  * 

lX#  4 

"zz  ' 

^  yy  )  cl 1  • 

M  =  Iyyq  ♦ 

iyyH  * 

^xx 

N  •  l22r  * 

lzir  * 

(lyy  ' 

rxx)P4- 

Ttiese  arc  the  equations  whicJi  describe  the  rotations. 

THE  FORCES  ANOjffljCjTC  ACTING  ON  THE  ROCKET 
HIE  AERODYNAMIC  FORCES  AND  MOMENTS 

Since  the  aerodynamic  forces  and  moments  depend  upon  the  missile's 
velocity  with  respect  to  the  surrounding  air,  preliminary  expressions 
must  be  obtained  before  these  forces  and  moments  can  be  computed. 

Let  Wy , ,  Wyi,  and  Nj,  be  the  components  of  the  wind  in  the  X'.Y’.Z* 
system.  The  wind  components  in  the  x,y ,z  system,  Wx,  W  ,  and  Wz,.arc 
obtained  by  applying  the  transformation  TXi2x.  This  can  be  conveniently 
written  in  matrix  form  as 

"w  1  fl.  1,1,1  f -cosait  sin*.  sinwt  -cos*.  sinwt  n 

X  1  2  3  j  L0  E0  * 

W  *  m  m  m  -sinait  -sin*L  cosut  cos*,  cosut  Wy,  (37) 

y  k  »  ■  !  o 

-WzJ  L  njn2n3  -1  L  0  C0S^Lo  iin*la  j  L  WZ  j  • 


The  velocity  components  of  the  missile  relative  to  the  wind  arc 
u'  ■  u-Wx ,  v’  ■  v-Wy,  w'  »  w-Wj,  and  the  relative  speed  is 


Va  -  l(u')2  ♦  (V)2  ♦  (w')2)  /2. 


There  are  several  angles  which  arc  used  to  calculate  the  aerodynamic 
forces  and  moments.  These  angles  are  shown  in  Figure  4.  They  c&n  be 
expressed  in  terns  of  velocity  components  as  follows;  the  angle  of 
attack,  a,  the  angle  of  sideslip,  6,  and  the  absolute  angle  of  attack,  6, 
given  by 


a  *  tan  • 


FI  GURU  4 


It  is  assumed  that  speed  of  sound,  V's,  and  atmospheric  density,  o, 
are  known  as  functions  of  altitude.  iow  the  Mach  amber  and  dynamic 
pressure  can  be  computed  from 


M.N.  -  JL  , 
Vs 


C43) 


and 


4' 


(44) 


Since  it  lias  beei.  assorted  that  the  ncrobyn.iric  forces  arul  r  orients 
can  be  cot ai neb  by  the  use  of  linear  aerodynamics,  the  forces  anil  nor.cnts 
a  re  [  S  ,  o  J 
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'y  "  l 


si  118  •  sina*]ii'h 

1  ‘  *’a  “'a 

P‘1 


1‘  a  I  -C,  sina*  *  t:;.  C™-)  sin&jq’S 

*•  '’u  I* 


L  -  [  C.  ♦  t  sin6] q'SU 

P  *  v  a 

M  -  [C  sina*  ♦  C  (3J-]  ♦  l_  (£")  sin8]  q’SU  , 


'raq  ‘2Va 


(45j 


•'  -  [‘Si!  sin8  *  C™  $->  ♦  S*  sin«*]q’Sd, 

a  H  a  “a  a 

where  S  is  the  reference  area  and  d  is  the  reference  diameter. 

nit  nirajsT  iorcl  and  nit  JtT  damping  moment 

Let  m  be  the  mass  of  the  rocket  including  the  unspent  fu*  1  an'1  let 
Am  be  the  change  in  mass  (due  to  burning  of  fuel)  during  a  small  time 
interval  At.  By  the  law  of  conservation  of  total  Momentum  we  can  equate 
the  momentum  at  a  time  t  to  that  at  a  time  t  ♦  At.  This  gives  us 

mv  •  (m  ♦  Am)  (V  ♦  AV)  ♦  Aa( V#  -  V)  (46) 

where  Vg  is  the  velocity,  relative  to  the  rocket,  of  the  exit  gases. 

By  dividing  both  sides  of  the  above  equation  by  At  and  then  taking 
the  limit  as  At  approaches  zero  we  get 

eV  *  -  m\'e,  (47) 

This  is  the  force  on  the  rocket  due  to  the  changing  moment  cm. 

Besides  -mVe  tiicrc  is  an  additional  force  due  to  the  differential 
in  the  pressure  at  the  exit  nozzle  and  the  atmospheric  pressure.  If  the 
pressure  at  the  exit  nozzle  is  Pe,  the  atmospheric  pressure  Pa,  and  the 
area  of  the  exit  nozzle  Ag,  then  this  additional  force  is  AgfPg  -  Py)  , 
giving  a  total  thrust  of 


T  -  mVe  ♦  (Pc  Pa)  Ac 


(4S) 


il  the  thrust  is  measured  .it  .1  test  stanu  (i.here  »  oj  in  ;m 
atmosphere  l*s  t  it  would  Iv 


Vt.  =  “VC  *  AC  C'c 

S’s  t  > 

(l-J) 

(  IS)  and  (-10)  we  get 

T  a  Ts.t,  *  Ac  (Pst 

■  l'a). 

(SO) 

'•e  have  assumed  that  the  thrust  acts  parallel  to  the  longitudinal 
axis  of  the  rocket,  hence  the  only  component  due  to  the  thrust  is  along 
the  x-axis  and  it  produces  no  moments. 

Since  a  rocket  rotates  about  a  transverse  axis  during  burning,  the 
gases  must  be  accelerated  laterally  as  they  flow  down  the  motor  tube. 

This  lateral  acceleration  produces  the  so-called  jet  damping  moment. 

The  following  expression  for  the  moment  is  derived  in  reference  [7] : 

Mp  *  m(lj  -  lfj)Wp, 

where  m  is  the  mass  flow  rate,  W  thn  instantaneous  pitching  velocity, 

1;  is  the  distance  between  the  vehicle's  Cg  and  the  exit  nozzle  and  lp  is 
the  distance  between  the  Cg  of  the  vehicle  and  the  propellant  Cg. 

The  components  of  the  moment  are,  for  a  symmetric  rocket, 

Mm  •  m(lj  -  l|)q 

Mj.  -  »(lj  *  l£)r  .  (52) 


nit  FORCE  UJC  TO  GIIAVITY 


Let  gs  be  the  average  value  of  g  at  sea  level  along  the  trajectory 
Then  the  value  of  g  at  an  altitude  h  is,  by  the  inverse-square  law. 


g  * 


(53) 


where  Uq  is  the  radius  of  the  earth.  Since  gravity  is  a  central  force 
it  acts  along  the  radius  vector  *  from  the  center  of  gravity  of  the 
missile  to  the  origin  of  the  X,Y,Z  system  and  hence  has  direction  cosines, 
in  the  X,Y,Z  system 


7rr~)'  (nr^r-ijf  (7 
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X2  ♦  YZ  ♦  Z2 


i>y  aprlyini:  the  Uaii.'h.watioi:  [\2x 


uc  a a vc 


fix  -  U.x  ♦  12V  *  1,:) 

•'  /X2  *  Y-  ♦  Z' 


I'  ‘  -/=^r=----~;~_(n.X  ♦  n  Y  *  n.  Z) 
>  /x2  *  Y2  ♦  :2  ‘  3 


s,  = 


i=_(n,X  ♦  n„Y  ♦  n.Z) 


(54) 


/72  ♦  Y2  *  Z2 

THli  EQUATIONS  OF  MOTION 

Tlic  six  equations  of  laotion  were  developed  in  the  preceding  section, 
and  the  forces  and  moments  which  were  assiT.ied  to  affect  tlie  motion  were 
discussed.  At  this  point  we  collect  the  previous  developments  so  that 
the  system  of  equations  may  be  presented  together  in  final  form.  The 
translational  equations  of  motion  are,  by  (31),  (45),  (49),  and  (54), 


mu  *  mrv  •  mqw  ♦  C»  i|  S  ♦  T  ■  nj;. 


(55) 


bi v  *  npw  •  aru  *  i  (--c  sin8 •  (-•;  j  (^r-)sina*  )  q'S  -  t.ify,  (5b) 
'  a  ”  a 


cim  *  mqu  -  rapv  ♦  {  Og  sina**C.-  (^r-Jsinfl  }  q’S  -  mgz 

4  ‘  «  «1 


(57) 


Ihe  rotational  equations  of  motion  arc,  by  (35),  (45),  (51),  and  (52), 
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I.  esc  equations  arc  mi.'-ruall;  uitij'iitc!  te  obtain  the  nut  ton  of  the 
rochet  in  the  x.y. z  system  At  the  sine  tire  liquations  ( 11.J  ,  (Is)'1,  ami 
(20)  arc  numerical  1/  ;nttgr atcii  to  iii>ra>!>  ’ ; ic  transformation  matrix  re¬ 
quired  to  express  the  results  in  the  other  systems.. 


SiJ'1 1AIIY 

A  six- degree -of- freedom  digital  siuulation  model  has  been  developed 
for  determining  the  trajectory  and  wind  effect  on  a  multistage  unguided 
rocket.  Although  it  was  not  specifically  mentioned  in  the  development 
of  the  equations  of  motion,  it  should  be  intuitively  obvious  tiiat  there 
is  a  discontinuity  in  these  equations  at  the  separation  and  expulsion 
of  any  booster. 

Several  assumptions  wore  made  m  the  beginning  of  tn*  development. 
Several  of  the  assumptions  could  be  easily  discarded  if  such  is  desirable. 
The  possibility  of  dropping  these  assumptions  will  be  discussed  below 

Assumption  No ■  1 .  Die  rocket  is  assumed  to  be  a  rigid  body  with  six 
degrees  ofTrcedoin.  It  appears  that  the  problem  would  be  unnecessarily 
complicated  if  the  rocket  were  not  assumed  to  be  a  rigid  body. 

Assumption  No.  2.  The  body  axes  are  principal  axes  of  inertia.  Inspee- 
tion  of  Equations  (23)  and  (34)  shows  that  this  assumption  is  clearly  not 
necessary.  It  was  made  for  two  reasons  Most  rockets  are  nearly 

symmetric  for  stability  reasons,  and  (2>  the  products  of  inertia  required 
in  lieu  of  this  assumption  are  not  usually  available.. 

Assumption  No.  3.  Linear  aerodynamics  arc  adequate  for  determining  the 
aerodynamic  forces  and  moments  Inspection  of  liquation  (45)  shows  that 
the  development  is  easily  adaptable  to  a  change  in  the  type  of  aerodynamics 
required. 

Assumption  No,  4.  The  earth  is  issuned  to  be  a  sphere.  This  assumption 
is  not  really  used  in  the  development  of  the  equations.  It  is  used  only 
in  the  analysis  of  the  results  of  .lie  computations.  It  would  be  simple 
enough  to  consider  the  earth  as  an  oblate  spheroid. 

Assumption  No.  5.  Gravity  follows  an  inverse  square  law.  This  assunp- 
tion  is  standarJ  for  the  types  of  missiles  considered. 

Assumption  No.  6  Hie  thrust  vector  acts  parallel  to  the  longitudinal 
axis  of  the  "rocket,  .‘lost  rockets  arc  d  signed  for  such  a  thrust  orien¬ 
tation.  If  it  were  desirable  to  consider  thrust  misalignments  for 
dispersion  analysis  it  would  be  a  simple  matter  to  compute  the  x,y,z 
component  of  the  tim.st  and  the  resulting  moments-  However,  for  appli¬ 
cations  indicated  m  this  repert  such  considerations  are  not  appropriate. 
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